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Thermal hysteresis in a micron-size Superconducting Quantum Interference Device (/x-SQUID), 
with weak links as Josephson junctions, is an obstacle for improving its performance for magne- 
tometery. Following the "hot-spot" model of Skocpol et al. [J. Appl. Phys. 45, 4054 (1974)] and 
by incorporating the temperature dependence of thermal conductivity of superconductor using a 
linear approximation, we find a much better agreement with the observed temperature dependence 
of the retrapping current in short superconducting Nb-based weak links and /i-SQUIDs. In addition, 
using the temperature dependence of the critical current, we find that above a certain temperature 
hysteresis disappears. We analyze the current-voltage characteristics and the weak link temperature 
variation in both the hysteretic and non-hysteretic regimes. We also discuss the effect of the weak 
link geometry in order to widen the temperature range of hysteresis-free operation. 



I. INTRODUCTION 

A micron-size superconducting quantum interference 
device (^-SQUID) consists of two superconducting 
Dayem bridges or weak links (WL) of dimension of 
the order of the superconducting coherence length, in 
parallel, forming a loop with area in the /zm 2 range. A 
single WL behaves very much like a Josephson Junc- 
tion [l| with the supercurrent approximately given by 
I = I c sin 9, where I c is the critical current and 9 is 
the phase difference across the junction. When two such 
junctions are fabricated in parallel in a SQUID, inter- 
ference between the two current branches gives an os- 
cillatory behavior of the critical current I c with the ex- 
ternal magnetic field @. The flux period is equal to 
the flux quanta $ = 2.05 XlO -15 T.m 2 . This makes 
the SQUID a very sensitive device to measure magnetic 
field. While the flux sensitivity achieved by conventional 
SQUIDs is better than 10" 7 $ /Vffi, for a /x-SQUID it 
has only been about KT 4 $ /VHz [|. /x-SQUIDs have 
been used to study the magnetization reversal of an 
isolated magnetic nano particle, the persistent current in 
phase- coherent rings [ij and also in scanning SQUID mi- 
croscopy Q . An improved sensitivity of /i-SQUIDs would 
be useful for probing ferromagnetic particles of smaller 
size or where the surface spins play an important role 

Other than the sensitivity, the hysteresis in /i-SQUIDs 
current-voltage (I-V) characteristic (see e.g. Ref. is a 
major hurdle as it (1) increases the measurement time, 
(2) complicates the measurement electronics, (3) changes 
the temperature of the sample placed in close proxim- 
ity with the /z-SQUID. Thus it is important to under- 
stand this hysteresis and devise ways of eliminating it. 
Hysteresis in the current- voltage characteristic is a very 
common phenomena for many superconducting nano- 
structured systems, especially WLs. It includes conven- 
tional Superconductor-Insulator-Superconductor (S-I-S) 
Josephson junctions 2], Superconductor-Normal metal- 
Superconductor (S-N-S) junctions 0, superconduct- 
ing nano- wires [§] and superconducting /x-bridges [il)l - [l2j . 



When the current is ramped up from zero across such 
junctions, the system suddenly switches to a non-zero 
voltage state at the critical current J c . After switching, 
when the current is ramped down, the system comes back 
to a zero-voltage state at a particular current, called the 
retrapping current I r . At very low temperature, the re- 
trapping current can be smaller than the critical current: 
I r < I c . This defines an hysteretic I-V curve. 

A number of models have been proposed in the last 
few decades to understand the hysteresis in supercon- 
ducting WLs. The resistively and capacitively shunted 
junction (RCSJ) Q model predicts the I-V curve for a 
conventional S-I-S junction very well. In this case, the 
capacitance across the junction is responsible for the hys- 
teresis. But for lateral junctions (either S-N-S junctions 
or constrictions) , the geometrical capacitance is too small 
to explain hysteresis. Hence an alternative theory of 
an effective capacitance C e tf was proposed (l2l |. where 
one equates the charge relaxation time R n C e ff with the 
Cooper pair relaxation time h/A. Here R n is the normal 
resistance and A is the superconducting gap parameter. 
The same method was recently extended to S-N-S junc- 
tions 0] by equating R n C e ff with the diffusion time of 
Andreev pairs. Though these methods reproduce some 
of the features of the I-V curves, no justification behind 
the origin of an effective capacitance has been found. 

Recently, Courtois et al. [13| have unambiguously 



shown, by directly measuring the electronic temperature, 
that heating is responsible for hysteresis in S-N-S junc- 
tions. According to the "hot-spot" model of Skocpol et 
al. [ll[, the heat generated in the resistive region of 
the WL raises locally its temperature above the criti- 
cal temperature T c . The temperature goes down to the 
bath temperature as one moves away from the hot spot. 
This gives rise to a normal metal-superconductor inter- 
face along the surface defined by T = T c . The interface 
location is self-consistently determined by the heat gener- 
ated and the coupling to the thermal bath. It was found 
that below a certain current, identified as the retrapping 
current, this interface becomes unsustainable and the WL 
turns fully superconducting. For a short WL, the hot 



spot may spread beyond the WL and into the electrodes. 
This "hot-spot" model reproduced most of the features 
of the I-V charac teristics o f superconducting WLs. It 
also predicted a y/1 — T/T c dependence of I r on T; how- 
ever the latter was not experimentally verified. Further, 
this model ignored the temperature dependence of the 
thermal conductivity of superconductor. Incorporating 
an approximate form for this temperature dependence, 
Tinkham et al. described the I-V characteristics of free 
standing superconducting nano- wires [3]. In this case, 
the N-S interface occurs inside the long nano-wire, mak- 
ing the problem one-dimensional. Again, this work did 
not include the temperature dependence of the hysteresis 
in the I-V characteristics. 

In this paper, we describe an effective one-dimensional 
thermal model to find out the temperature profile near a 
short WL connected to wide electrodes. We calculate the 
I-V characteristics as well as the (bath) temperature de- 
pendence of the retrapping current. Our model predicts 
how the normal-superconducting (N-S) interface position 
varies with various parameters like temperature, current, 
and geometry. We also discuss the detailed temperature 
profile and how it changes with the bias current. Us- 
ing the temperature dependence of the critical current 
near T c , we find that above a certain temperature Th, 
hysteresis disappears. The effect of the WL geometri- 
cal parameters on Th is discussed quantitatively. As the 
same model is directly relevant to /i-SQUIDs, we test it 
on several such samples. Our model fits our data very 
well. Finally, we discuss how the non-hysteretic regime 
can be achieved over a wider temperature range, followed 
by conclusions. 



II. THERMAL MODEL OF HYSTERESIS 

Following the "hot-spot" model [Tl| , we consider a sin- 
gle WL connected to two extended electrodes, as shown 
in Fig. [T^,, and investigate the temperature distribution 
around it in the resistive state. We assume a local quasi- 
equilibrium condition so that a local temperature can 
be defined at each point of the sample. The length and 
width of the WLs under study are in 50-200 nm range, 
i.e. of the order of the coherence length of bulk Nb (£jV6 ~ 
39 nm). In this range, a WL behaves very much like a 
Josephson Junction 1]. Since the WL size is very small, 
we assume that in the resistive state the entire WL re- 
gion stays at a uniform temperature. In reality, the WL 
will have certain spatial temperature variation, but what 
will matter here is the heat evacuated out of the WL. 
This assumption will not affect our conclusions as long 
as the WL temperature is not so large as to affect its 
resistance. We also assume the pads connecting to the 
WLs to be much wider than the length scale over which 
the temperature relaxes to the bath temperature. 

At a given bias current, we can divide the device into 
three regions as shown in Fig. [Ud: 1) the narrow WL at 
a uniform temperature T\ consisting of the rectangular 
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FIG. 1: (a) Scanning electron micrograph image of a WL. 
(b) Sketch of the sample geometry with the three regions dis- 
cussed in the text. 



(width w and length I) area in the center and terminating 
into a semicircle of radius r\ = w/2 at each end, 2) a nor- 
mal state (T > T c ) semicircular region in each electrode 
with ro > r > w/2, and 3) a superconducting (T < T c ) 
region for r > tq. We have assumed a rounded WL of 
radius r\ = w/2 to avoid logarithmic divergence in the re- 
sistance calculation. This approach is unlike the hot spot 
model for a long WL [ll|, where the hot-spot develops 
near the center of the WL. Beyond the WL, we assume 
that the heat conducts away radially in the bulk of the 
film. Thus the temperature also decreases radially inside 
the two electrodes, reaching the bath temperature T b far 
away from the WL. This enables us to use an effective 
one-dimensional model for finding the spatial variation 
of temperature. The source of heat is the resistive dissi- 
pation in the normal region, which extends up to a radial 
distance ro in each electrode, thus defining a N-S inter- 
face with T = T c between normal and superconducting 
regions. 

In the resistive state, the Joule heat near the WL re- 
gion is removed in two ways: conduction within each elec- 
trode and surface heat flow from film's bottom surface to 
the substrate. The latter is assumed to be proportional 
to the temperature difference between the film and the 
substrate. This approximation has been used extensively 
[Til [TBI [l6|: we will discuss its validity later. We neglect 
the heat loss from the top surface of the film as we oper- 
ate in a vacuum cryostat. We also assume that the entire 
substrate stays at the bath temperature T b - Thus, the 
general heat flow equation can be written as: 

dT r 

-n—2Trrd+ a(T - T b )2nr dr + a(T x - T b )A 
dr J ri 

= I 2 R{r). (1) 

Here k, is the thermal conductivity, d is the thickness 
of the film, A is the area of the WL region given by 
A = (wl + irri 2 ), I is the current. The surface heat loss 
coefficient a is expressed in W/m 2 K. R(r) is the resis- 
tance of the electrode within a radius r including the re- 
sistance of the WL. Assuming a radial and isotropic cur- 
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rent flow in the electrodes giving rise to circular equipo- 
tential lines, the resistance R(r) is given by: 



R(r) = R + R c \n(- 



(2) 



where nR c /2 = pn / d is the square resistance of the film, 
with p^r as the normal state resistivity. Thus R c defines 
a characteristic resistance of the film. The resistance Rq 
of the WL is given by Ro = (I + 2ri)p^ /wd. 

Since we consider low temperatures, we ignore the 
phonon's contribution to the thermal conductivity. In 
the normal state, the electronic part of the thermal con- 
ductivity be found using the Wicdemann-Franz 
law: K n = LqT/pn, where K n is the thermal conductivity 
in the normal state and Lq is the Lorentz number. In the 
superconducting state, as the number of quasi-particles 
decreases significantly, one expects the thermal conduc- 
tivity n s to be much less. At very low temperature, since 
only very few quasi-particles are left to carry thermal en- 
ergy, k s can be exponentially small. We use here a linear 
approximation, k s /k„ = T/T c , which gives k s = n n at 
T = T c as expected. From the theory [I?), this linear 
approximation is well justified near T c . The same ap- 
proximation was also used by Tinkham et al. [14| . 

Using the above expressions for k and R(r) and differ- 
entiating Eq. [JJwith respect to r, one gets: 
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The boundary conditions are: 1) at r = ri, T 
Eq. m gives - K (T 1 )^2^r 1 d + a(Tj - T b )A = 
at r = ro, T = T c , T and ^ are continuous and 3) for 
r — > oo, T = T b . The radius r and WL temperature 7i 
have to be found self-consistently using these boundary 
conditions. 

An inspection of the above two equations gives us a 
length scale, 



and a current scale 
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Here Io would determine the scale of the retrapping cur- 
rent I r , while rj would determine the length scale of tem- 
perature variation. For WLs based on a Nb film de- 
posited on a Si substrate, one typically uses a thickness 
of 20 to 150 nm. Depending upon the detailed prepa- 
ration method, some typical parameters would be Pn_ = 
15-50 pfl cm, T c = 6-9 K and a = 1-3 W/cm 2 K [lj]. 



Using L Q = 2.44xl0~ 8 W.il/K 2 , we get t] -1-3 pm and 
I - 0.5-2 mA. 

Eq. [3] and [4] can be written in terms of the dimension- 
less variables x= r/rj, t = T/T c , i = I/Iq, x\ = fi/r], xq 
= r^/rj and tf, — T(,/T c as follows: 
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- (t - t b ) = (xi <x< x ), (7) 
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dx 



(t-t b ) = (x>x Q ). (8) 



In terms of the reduced variables, the boundary condi- 
tions become: 1) at x = Xi, -xihjt + (*i ~ = 

TtdRg '2 



2p N 



2) at x = Xq, t = 1 and dt/dx is continuous, 
and 3) for x — > oo, t = tf,. For short weak links, us- 
ing 2^rr << 1, the first boundary condition becomes 

*ii = -0**/*i with 13 = ^ = jg = f(l + hi 

Eq. [7] and [8] are second order and non-linear differen- 
tial equations that can be solved only numerically. We 
wish to go beyond the approximation of k being indepen- 
dent of temperature, which would give solutions in terms 
of modified Bessel functions as discussed by Skocpol et 
al.[ll(. We choose to simplify the above equations by 
substituting yi = t 2 and y2 = t 3 in Eq. [7] and respec- 
tively, yi and ?/2 then satisfy: 
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1- 2{y/yl-t b )x (xi < x < x ), (9) 
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= 3(y 2 / - t b )x (x > XoX-lO) 



Let us first consider the superconducting region (x > 
xo) described by Eq. [TO] In this equation, j/2 varies be- 
tween t 3 and 1 . For this range of j/2 , we linearly approx- 

1 /3 

imate the (y 2 — H) term as: 
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(11) 



so as to keep the end points of (y 2 — h), i.e. at t = t b 
and (1 — t b ) at t = 1, fixed. This approximation be- 
comes more and more accurate as t b approaches 1, i.e. 
the bath temperature T b approaches the critical temper- 
ature T c . Eq. [TU] then reduces to the modified Bessel 
equation s 2 £" + s£' — s 2 £ = 0, where £ = y 2 — t 3 and 
s = Xx with A = ^3/(1 + t b + i 2 ). With the boundary 
condition t = t b (i.e. £ = 0) at x — > oo, the only accept- 
able solution is £ = CK (s), where K is the modified 
Bessel function of second kind and zero degree. Using 
the boundary condition t = 1 at x = Xq, we get the final 
solution for x > Xg as: 



t' 



1 



t 3 



K (\x 



K (Xx). 



(12) 



Fig. [5] shows for comparison the numerical solution of 
the non-linear Eq. [TU] and the corresponding solution to 
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FIG. 2: Comparison between exact (Eq. 1111 black) and ap- 
proximate (Eq. 1121 red) solution for the reduced temperature 
profile for x > xq. The respective calculated slopes at the 
origin are - 5.1 and - 4.80. Parameters are tb = 0.5 and 
x = 0.015. 



the linearized Eq. [T^fbr tb — 0.5 and xq — 0.015. The 
overall shapes of the curves are similar, justifying our 
approximation. 

Let us now consider the normal region (x < xq). Eq. 
is difficult to linearize as yi varies between t b and t\, and 
t\ is not known before-hand. The nature of boundary 
conditions does not allow a simple numerical solution. 
We make the approximation to neglect the surface loss 
term, i.e. the (y/yx — i&) term in Eq. [HI This is justified 
for finding the retrapping current I r in the regime xo > 
X\, in which case the heat loss to the substrate is not 
significant as compared to the heat conducted out. With 
this approximation, there is an analytical solution: y\ — 
—i 2 [(In a; + Ci) 2 ) + Ca] • Here C\ and C2 are constants 
to be found from the boundary conditions: y% = t 2 = 1 



at x = xo and dyildx = 

X\ < X < X : 



-2/3i 2 jx\. Finally, we get for 



t 2 = 1 _ e [(In - + pf - (In ^ + P) 2 ] . (13) 

X\ X\ 

This relation gives the temperature profile for x < xq 
and determines the WL temperature t-y = t(xi) in terms 
of xq. To find Xq, we have to use the continuity of dt/dx 
at xq using solutions given by Eq. [12] and [13] This gives 
the following transcendental equation for xq: 



\xq{1 ~ tpK^Xxo) 
Z[\n{x a /x 1 ) + 0\K o (\x o ) 



(14) 



As shown in Fig. [3Ji, the right hand side of above Eq. 
1141 features a minima in current i as a function of xq. 
This means that below this current, the Joule heat is not 
sufficient to uphold a normal metal-superconductor (N- 
S) interface with T = T c . This current is thus identified 
as the retrapping current i r = I r /Io- Fig. 13b shows that 
it decreases with increasing bath temperature, whereas 
the related xq increases. At high temperature, a regime 
where the retrapping current exceeds the critical current 
(i r > i c ) can be reached. In this case, the WL is resistive 
while its temperature stays below T c . Only if the bias 
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FIG. 3: (a) Plot of Eq . 1 141 right-hand side as a function of xo 
at different bath temperatures. The minima shown by arrows 
define the retrapping current i r . (b) Variation of the retrap- 
ping current as a function of the bath temperature, (c) I-V 
characteristics near i r at three different bath temperatures as 
indicated in the figure. The current and the voltage axes are 
normalized with respect to Io and IoR c , respectively. All the 
curves are plotted for /3 = 3.5 and x\ — 0.015. 



current becomes larger than i r , does a N-S interface with 
T = T c appear at tq (> ri). We will discuss this point 
in more detail later. 

Eq. [T3] provides the relation between the current bias / 
and the N-S interface position xq. One can then calculate 
the resistance R(xq) using Eq. [2] The related current- 
voltage characteristic V = IR(xq) is plotted in Fig.[5J: for 
different bath temperatures. At low voltage, a negative 
differential resistance branch appears, meaning that, in 
this regime, for a given current, the voltage can have 
two distinct values. Since this branch is unstable under 
current biasing 18], only the higher voltage is accessible. 
But if one performs voltage-biased measurements, then 



5 



one can access the negative differential resistance branch 
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m micro- 



as well, as was observed by Skocpol et al. 
bridges and Steinbach et al. [l^ in Josephson junction. 

For illustration, let us now consider a WL biased at 
its retrapping current. Fig. |4] shows the radial tempera- 
ture distribution for different bath temperatures and for 
some typical values of (3 and x\. Expectedly, at large dis- 
tance, the temperature asymptotically decreases to the 
bath temperature. The temperature profile close to the 
WL exhibits a large temperature gradient as compared 
to the superconducting region, see Fig. |4j;, d. The inter- 
cepts of the different curves with the dotted horizontal 
lines representing t — 1 indicate the location of the N-S 
interfaces. The temperature values at x = x\ (= 0.015 
here) indicate the WL temperature. 

Still at retrapping, Fig.[5^,b shows the variation of the 
N-S interface position ro in units of r\ as a function of the 
bath temperature for different X\ and /3 values and as a 
function of ft at a fixed bath temperature for different x\ 
values. For large values of (3, i.e. for long WLs, the N-S 
interface is closer to the WL. Fig.JSJ; and d show the tem- 
perature of the WL as a function of the bath temperature 
for a combination of f3 and x\ values. We observe a non- 
monotonic behavior, which is due to the increase of the 
thermal conductivity with increasing temperature. The 
contrast between, on one hand, the monotonic evolution 
of the current i r and the N-S interface position x$ with tj, 
and, on the other hand, the non-monotonic evolution of 
the WL temperature at retrapping indicates that it is the 
size of the normal region and not its local temperature 
that governs the amplitude of the retrapping current. 



III. TRANSPORT EXPERIMENTS ON WEAK 
LINK /^-SQUIDS 

We have tested the above model on /i-SQUID samples. 
A micrograph of one such device is shown in Fig. |6]inset. 
We use Nb films deposited using DC magnetron sputter- 
ing in a chamber with a base pressure in the 10 -7 mbar 
range. For most of the samples, Nb thin films are de- 
posited on a Si wafer and a photo resist is spun on the 
films. Using optical lithography, we form a coarse pat- 
tern (several /im size) on this resist, which is transferred 
to the film by wet chemical etching using dilute hydroflu- 
oric acid (HF). The final desired pattern is obtained by 
finer milling with the help of Focused Ion Beam (FIB). 

The film thickness d was measured using a profilometer 
across a step made by masking during deposition. The 
width w and length £ of the WLs were estimated from 
the SEM images. The sample thickness varies between 
30 to 75 nm, whereas the width and length of the WLs 
vary from 50 to 200 nm. All the WL dimensions are thus 
much smaller than the length rj defined earlier. For most 
of our devices, the maximum asymmetry between the two 
junction is less than 10% both in length and width. 

In this article, we report on four samples whose de- 
tailed parameters are given in Table HI Transport exper- 
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FIG. 4: Temperature (normalized) evolution with radial dis- 
tance (normalized) in a WL biased at its retrapping current 
at different bath temperatures indicated in the figures. The 
parameters are x\ — 0.015 and (a) (3 — 3.5, (b) /3 = 1.5. 
(c) and (d) are magnifications near the normal region, corre- 
sponding to (a) and (b) respectively. The intersection with 
the dotted lines at t = 1 indicate the N-S interface position. 
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FIG. 5: Variation of the N-S interface position ro in units 
of n at the retrapping current as a function of: (a) bath 
temperature if, for different values of x\ and f3, (b) parameter 
j3, at tb = 0.15 for different values of xi. (c) and (d) show the 
variation of the WL temperature at retrapping as a function 
of bath temperature for (c) /3 = 3.5 and (d) j3 — 1.5 for xi = 
0.015 and 0.02. The arrows show the maxima. 
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FIG. 6: Resistance vs temperature curve for Sample 2 down 
to 4.2 K at a bias current of 0.1 mA. Inset shows the SEM 
image of a typical /^-SQUID with a loop area 3.5 x 3.5 /im 2 . 



iments were performed down to 300 mK in a 3 He cryo- 
stat. We used r.f. niters at several stages of the cryostat 
to minimize noise. Measurements were done in current 
bias mode using a d.c. current source. No magnetic field 
was applied for the data presented here. 

Fig. |6] shows the temperature variation of Sample 2 
resistance down to 4.2 K. The main and sharp transi- 
tion with an onset at 7.44 K is expectedly for the bulk 
film. The other transitions (steps) correspond to rela- 
tively smaller pads connected to the WLs. Since the re- 
sistance has a large tail, it is difficult to define the critical 
temperature T c from this data. We therefore define T c 
from I-V measurements (discussed below) as the temper- 
ature above which I c is zero, i.e. the I-V curve is ohmic. 

Sample 1 I-V characteristics at different bath temper- 
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FIG. 7: (a) I-V curve of Sample 1 at four different bath tem- 
peratures. The plots of 1.5 K, 4.5 K and 5.5 K have been 
shifted upwards by 5, 10 and 15 mV, respectively, for clarity, 
(b) Temperature dependence of the critical and retrapping 
currents for the same sample. 



atures are shown in Fig. [7^. At low temperature, the 
curves are clearly hysteretic. From this type of data, we 
experimentally define the critical current as the maxi- 
mum current up to which no measurable voltage is ob- 
served when the current is ramped from zero. Here, we 
do not distinguish between critical current and switching 
current. In the retrapping branch, we define the retrap- 
ping current as the current at which the resistance goes 
back to zero. For most of the samples, with the above 
definitions, the detection of I c and I r are accurate within 
about 1% for T b < 1 K and about 10% near T c /2. Close 
to T c , because the transition region width, I c or I r can- 
not be defined with an accuracy better than 50%. Fig. 
[7Jd shows the variation of critical I c and retrapping I r 
currents as a function of bath temperature for the same 
sample. Above a temperature Th, the retrapping and the 
critical currents are equal, meaning that hysteresis in the 
I-V curve has disappeared. 

Fig. [8^ and b show the experimental I —V at the low- 
est temperature for the two samples together with their 
fit by our model. The current and voltage are normal- 
ized with respect to the fit-derived parameters Iq and 
R c Iq- In Fig. [SJ: and d, the variation of retrapping cur- 
rent as a function of temperature are shown. Here also 
the current is normalized with respect to Iq, whereas the 
temperature is normalized with respect to the critical 
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TABLE I: Comparison between various experimental (exp) and fit parameters: WL length £ and width w, film thickness d, 
characteristic resistance R c , critical temperature T c , current scale 7q. 



Sam. 


I 


w 


d 


Rc (exp) 


Rc (fit) 


T c (exp) 


T c (fit) 


J Q (exp) 


Io (fit) 


no. 


(nm) 


(nm) 


(nm) 


(n) 


(SI) 


(K) 


(K) 


(mA) 


(mA) 


1 


95 


75 


45 


2.7 ± 0.3 


2.5 ± 0.5 


5.50 ± 0.10 


5.25 ± 0.53 


0.64 ± 0.10 


0.81 ± 0.10 


2 


100 


100 


55 


2.3 ± 0.2 


1.6 ± 0.3 


5.70 ± 0.10 


5.60 ± 0.56 


0.77 ± 0.12 


0.78 ± 0.09 


3 


150 


145 


65 


2.0± 0.2 


1.9 ± 0.3 


5.80 ± 0.10 


5.60 ± 0.56 


0.91 ± 0.14 


0.97 ± 0.12 


4 


150 


150 


30 


5.3 ± 0.5 


5.1± 0.5 


4.5 ± 0.10 


4.35± 0.44 


0.27 ± 0.04 


0.24 ± 0.03 



temperature T c . The values of the fit parameters R c , Iq 
and T c together with the experimental parameters are 
listed in Table HI Here the experimental R c = 2p n /ird is 
calculated by measuring the resistance of a known rect- 
angular geometry and Iq — 2T c ^/Lq/R c is calculated by 
using the value of above R c and experimental T c . We 
have obtained a similar agreement with Sample 3 and 4 
(not shown here), whose experimental and fit parameters 
are also included in Table Q] We have also used the same 
model for single WLs and /x-SQUIDs from Hasselbach et 
al. Q with a different geometry. In both cases, we could 
fit both the low temperature hysteretic I-V curves and 
temperature dependence of the retrapping current with 
our model very well. 

Fig. [SJI also shows a fit of the retrapping current with 
the VI — t dependence from Skocpol et al. [llj . The 
fitted coefficient 0.18 compares reasonably with the esti- 
mated value of 0.34, assuming a Wiedemann-Franz law 
to get the normal-state thermal conductivity. Neverthe- 
less, our model gives a clearly much better agreement, 
which we attribute to incorporation of the superconduc- 
tor thermal conductance temperature dependence in our 
model. 



IV. DISCUSSION 

The exact thermal model for our system is quite in- 
volved with complicated non-linear differential equations. 
In this paper, we have tried to simplify them in a way 
that the essential features are preserved. This simplified 
model fits the experimental data very well. Nevertheless, 
several approximations need further discussion. 

We assumed the width of the connecting pads to be 
much greater than 77, but in actual experiments it is com- 
parable to it. Therefore the actual thermalization would 
be poorer than what is being assumed; we may be slightly 
overestimating x\. As it is difficult to estimate a for our 
samples and it actually has a temperature dependence 
[l5| . the determination of 77 and hence x\ can again be 
significantly erroneous. For most of the cases, we could 
fit our data with a 20% variation in the value of X\ by 
adjusting the other parameters. 

We have made the hypothesis that the electron and 
phonon temperatures are equal in the superconducting 
region. The electr on-p honon coupling power in a vol- 
ume V is given by go|, P = Y,V(Tj? - T p 5 ), where E = 




0.00 



0.0 0.2 0.4 0.6 0.8 1.0 



0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 8: Experimental (black) and numerical fit (red) of nor- 
malized I-V curves for Samples 1 and 2 at Tb — 300 mK with 
(a) P = 3.5 and x\ = 0.015 for Sample 1 (b) /3 = 3 and 
Xi — 0.02 for Sample 2. (c) and (d): Variation of the normal- 
ized retrapping current with the normalized bath temperature 
for Sample 1 and 2. The black dots are the data and the red 
curves are fits based on Eq. 1141 The fit parameters are listed 
in Table IT] The blue curve in (d) is fitting with the Skocpol 
et al. [lit (their Eq. 14) prediction i r = 0.18(1 - t) 1/2 . 



3 K 5 is the electron-phonon coupling pa- 
Taking typical experimental values I r 



2.4 x 10 9 W.nT 
rameter for Nb. 
= 0.12 mA, R c = 2.4 ft, r /n = 1.1 (that gives R{r ) 
= 8 ft), we get a dissipated power P = I^R(r ) = 0.12 
/j,W. Nearly all of this resistive heat is transmitted to the 
substrate in the superconducting region only. Though 
the temperature decreases sharply making the heat loss 
rather non-uniform, the effective size of this region is of 
order 77, which ranges between 1 and 3 /im. Taking 77 = 
2 /jm, d = 50 nm and an average electron temperature 
T e = 4K, the volume of the superconducting region is, 
V = Trrj 2 d = 3.5 x 10~ 19 77i 3 . This gives the temperature 
difference T e —T p as 0.06 K only, validates our hypothesis. 

We took a linear approximation for the surface loss 
term. The metal film and the substrate phonons ex- 
change heat through a Kapitza resistance, giving a power 
K a (T 4 — T£) per surface unit, with K a as the Kapitza 
constant [2(|. For T close to the above expression 
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can be approximated as 4if a T b 3 (T — Tf,). From the tem- 
perature profile in FigHl we can say that for most of 
superconducting region the above approximation is valid 
except for the region close to the N-S interface (r ~ ro). 
However, if the bath temperature is close to the critical 
temperature then for the entire WL the above approxi- 
mation would be valid. 

We also neglected the heat loss to the substrate from 
the normal region, i.e. the WL and the semicircular re- 
gion between 77 and r . Let us compare the heat transfer 
to the substrate P s and the heat conducted out P c un- 
der the linear approximation. Considering only one half 
of the film, we can approximately write P s = awr^ti — 
tb)T c /2. We can also write P c = — K„7rrod(^) r=ro = 
irrodLoT^ti — l)/(p„(ro — ri)). Here we assume a linear 
temperature decrease within the normal region, which 
is fairly justified according to Fig. Eh, d. For ti = 1.1, 
t b = 0.05, n = 100 nm, r = 2ri, a = 5 W/cm 2 .K, 
d = 50 nm, p n = 25 fifl cm, and T c = 8 K, one gets P c = 
0.2 fiW and P s /P c — 0.1, which confirms our assumption. 

The surface heat loss from the WL normal state region 
was neglected by assuming (ii — h)^^ "C /3i 2 . Let us 
check the argument for the worst possible case; i.e. at 
lowest possible temperature and for longer WLs. Taking 
t b = 0.05, t x = 1.1, A = 300 x 300 nm 2 , rj = 1 /jm at 
i = i r = 0.15, one gets (ti — h)^^ — 0.015, whereas 
with (3 = 7r (which corresponds to i — w — 300 nm) and 
i = i r = 0.15 we get f3i 2 = 0.072, which is almost 5 times 




3.0 3.5 4.0 4.5 5.0 



T b (K) 



FIG. 9: Variation of critical current with the bath tempera- 
ture in the high temperature regime for Sample 1. The blue 
line is a straight line fitting I c = I c q(1 — Tb/T c ) with I c q = 0.84 
mA and T c = 5.05 K. 



0.15 



C 0.10 
q 0.05 



0.00 



\ /?=3.5,x 7 = 



0.015 



X. 




higher than (ti — tb) 



0.75 0.80 0.85 0.90 0.95 1.00 
tu 



V. WHEN DOES HYSTERESIS DISAPPEAR? 

A key feature is the disappearance of hysteresis at high 
temperature. In general, there is a particular bath tem- 
perature Th at which I r and I c are equal. Above this 
temperature, for I > I c the current is large enough to 
kill the superconductivity in the WL, making it resistive. 
But the related Joule heating is not sufficient to raise 
the WL temperature above T c and provide an N-S in- 
terface. In order to find the crossover temperature Th, 
we need an expression for the temperature dependence 
of the critical current. For a bath temperature near the 
critical temperature, we can use the expression: 0, [2l| 
I c R n = ~fT c (l — tb), where 7 = 635 fiV/K and R n is the 
normal state resistance. In practice, 7 can vary signif- 
icantly. Taking R n — Rq, one can simplify the above 
equation to: 

i c = — t=(1-4). (15) 

Here, i c = I c /Io and we have used 7 = 635/iV/K, Lo — 
2.44 x 10~ 8 W.fJ/K 2 . In Fig. El we plot the variation of 
I c with the bath temperature for Sample 1 above T c /2. 
From the linear fit, we extract 7 = 930/iV/K with the 
above L value, (3 = 3.5 and R c = 2.2851. 

In Fig. [TU1 we plot the variation of i c (red curve) and 
i r (black curve) as a function of the normalized bath 



FIG. 10: Temperature dependence of i c (red) following Eg. 1151 
and of i r (black) calculated from our model, near the critical 
temperature. The parameters used for the plot are given in 
the figure. At tb = th, the two curves cross each other, so 
that hysteresis disappears for higher temperatures. 

temperature t b near tb = 1 using Eq. [15] for ft = 3.5 
and X\ = 0.015. The crossover temperature th is then 
straightforwardly determined from the intercept of the 
two curves. Let us point out here that the critical cur- 
rent I c and the retrapping current I r are controlled by 
two different physics, with I c dependent on the WL su- 
perconducting properties and I r on the heat dissipation. 
This justifies that these currents have a different temper- 
ature dependence. 

In Fig. [TTJ we plot the variation of th as a function of 
j3 for three different values of x\. The top axis refers to 
i/w, which is related to /3 by the formula: (3 = f (1+ — )• 
The upper limits (3 max of the parameter (3 are chosen in 
a way that at this point xq = 27, i.e. the N-S interface is 
at the WL boundary, as beyond this value our short WL 
approximation does not hold. From Fig. II H this occurs 
at a th of about 0.8. 

In the hysteresis- free regime t > th, the detailed tem- 
perature profile i r > i > i c can be found by solving 
Eq. [TO] for x > x\, i.e. the superconducting region. The 
boundary conditions used for solving this equation are 1) 



l/w 

0.9 1.2 1.5 1.8 2.1 2.4 



0.90n 

0.87 
^ 0.84 X 



xj= 0.005 



0.81 



0.78 



x,= 0.03 xj= 0.015 



3.0 3.6 4.2 4.8 5.4 

P 

FIG. 11: Variation of the hysteresis crossover temperature th 
with the parameter f3 for three different values of x\ (i.e. WL 
width). The top axis represents l/w, calculated using the 
formula: /?= f(l+ £). 



at x 



TrdRpi 
2pw 



and 2) t = tf, as x — > oo. 
Solutions for t were found numerically and are plotted in 
Fig. na. 

Fig.lPTb shows the temperature ii of the WL as a func- 
tion of current for two bath temperatures, tb — 0.90 and 
tb = 0.75, respectively above and below the hysteresis 
threshold th = 0.83. As the current is ramped up from 
zero, the WL temperature jumps from t\ ~ tb to a higher- 
value at i = i c - For a bath temperature above the hys- 
teresis threshold (tb > th), there is another jump in WL 
temperature at i = i r and after this the temperature 
keeps on increasing. The t\ vs i behavior remains same 
when i is ramped down, i.e. there is no hysteresis. Let 
us point out that there is no actual retrapping at this 
i r value, but the appearance or disappearance of a N-S 
interface close to the WL. For a bath temperature below 
the hysteresis threshold (tb < th), the behavior shows an 
upward jump from t\ = tb to a higher value when the 
current is ramped up through i = i c , and a downward 
jump to t\ = tb when current is ramped down through 
i — i r . Hysteresis is thus present. 

Fig. [T2b shows the i — v curves as calculated from the 
location of the N-S interface (above i r ) and the resistance 
of the WL (below i r ) at the same two bath temperatures 
above and below th ■ We see a close resemblance between 
the temperature and voltage curves as a function of bias 
current. While below th, the i — v spectra describe well 
the experimental curve, above th the calculated curve 
shows an extra step at i r arising from the sudden cre- 
ation of the N-S interface at a position ro > r\. This 
step was not observed in experiments. We believe that 
the predicted extra step may get significantly rounded 
as the WL temperature approaches T c . The supercon- 
ducting region outside r\ will be close to T c , reducing 
its critical current density. The exact shape of the I-V 
curve in the non-hysteretic regime will then be dictated 
by thermally activated phase slips 0, [23j for i < i c and 



0.96 



0.92 



0.88 



0.84 



\ /= 0.076 


(a) 






0.024 




0.00 0.04 0.08 0.12 0.16 

i 

1.2i 




0.00 0.04 0.08 0.12 0.16 
I 

FIG. 12: Calculation results with parameters /3 and xi equal 
to 3.5 and 0.015 respectively, which give th= 0.83. (a) Cal- 
culated temperature distribution at a bath temperature tb = 
0.85 above the threshold temperature th for normalized bias 
currents of 0.024, 0.054 and 0.076. Here i r is 0.093. (b,c) 
Variation of WL temperature (b) and i — v curve (c) as a func- 
tion of bias current for bath temperatures above the threshold 
tb = 0.75 (black) and below tb = 0.9 (red). The blue arrows 
indicate the direction of current sweep. 



superconducting fluctuations for i > i c [23J. 

Let us now consider how one could manipulate the 
crossover temperature th- At a fixed x\ (which is propor- 
tional to the width w), th decreases with the increase of (3 
(and hence the length £), see Fig. [TTJ This is desirable if 
we want hysteresis to disappear at low temperature. But 
the adjustment of j3 to any arbitrary value is impossible, 
since we wish the WL to behave like a Josephson Junc- 
tion, which implies the condition: w < I ~ £ This 
gives a lower bound on (3 and x\ and hence th in general. 
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Therefore, it is generally not possible to eliminate hys- 
teresis for these WL junctions just by manipulating the 
WL width and length. 

However, since j3 is the ratio of WL resistance i?o to 
the characteristic resistance R c of the film, we can effec- 
tively increase j3 by increasing the WL resistance. This 
can be done by reducing the thickness of the WL alone. 
The reduction of the whole film thickness (d) including 
the connecting electrodes can also reduce T^: in this case, 
•q oc y/d is smaller, giving a larger xq. This makes hys- 
teresis disappear at smaller temperatures (see Fig. [11]) . If 
the critical temperature T c is not affected, a wider tem- 
perature span for the hysteresis-free regime is achieved. 
This improvement was observed by Tinkham et al. [H[ 
for superconducting nano-wires. With an appropriate 
choice of substrate and growth conditions, the resistivity 
and critical temperature of very thin films can remain al- 
most unaffected by thickness reduction (for Nb see Ref. 
Hil), enabling similar results to be obtained with contin- 
uous films. Another possible way is to reduce the critical 
current I c (and possibly increase j3) by either making the 
WL with a poor superconductor or completely replacing 
it by a normal metal. Angers et al. [7| were able to get 
a T h < 1.2K for ^-SQUIDs made with SNS type weak 
links based on Nb. 



VI. CONCLUSIONS 

In conclusion, we have described a thermal model for 
understanding the hysteresis in the I-V curve of short 



superconducting WLs and their extension to /i-SQUIDs. 
Using this model, we have calculated the detailed I-V 
characteristics and the temperature profile near the WL 
as a function of bath temperature. We have obtained 
a good agreement between experiments and theory in 
terms of I-V characteristics and their temperature de- 
pendence. A key finding of this model, which again is in 
agreement with the experiments, is the disappearance of 
hysteresis above certain temperature. We have discussed 
how one can adjust the WL geometry in order to widen 
the temperature range of this hysteresis-free regime of 
WL-bascd ^-SQUIDs. 
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